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Abstract 

We calculate differential cross sections for exclusive production of heavy charged scalar, weakly 
interacting particles (charged Higgs bosons, charged technipions, etc.) via photon-photon ex¬ 
changes in the pp —> ppH H reaction with exact 2 —> 4 kinematics. We present distributions 
in rapidities, transverse momenta, and correlations in azimuthal angles between the protons and 
between the charged Higgs bosons. As an example, the integrated cross section for i/s = 14 TeV 
(LHC) is about 0.1 fb and about 0.9 fb at the Future Circular Collider (FCC) for n/s = 100 TeV 
when assuming m H ± = 150 GeV. The results are compared with results obtained within standard 
equivalent-photon approximation known from the literature. We discuss the role of the Dirac and 
Pauli electromagnetic form factors of the proton. We have also performed first calculations of 
cross sections for the exclusive diffractive Khoze-Martin-Ryskin mechanism. We have estimated 
limits on the ghH+H- coupling constant within two-Higgs dublet model based on recent experi¬ 
mental data from the LHC. The diffractive contribution is, however, much smaller than the 77 
one. The Z 7 , 7 Z, and ZZ exchanges give even smaller contributions. Absorption corrections 
are calculated for the first time differentially for various distributions. In general, they lead to a 
damping of the cross section. The damping depends on the M^+h- invariant mass and on t four- 
momentum transfers squared. In contrast to diffractive processes, the larger the collision energy, 
the smaller the effect of absorption. We discuss a possibility to measure the exclusive production 
of two charged Higgs bosons with the help of so-called "forward proton detectors" at the LHC 
experiments. 
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I. INTRODUCTION 


There are several reasons why exclusive reactions are interesting [jl, 2], One of them 
is the possibility to search for effects beyond the Standard Model (SM). The main advan¬ 
tage of exclusive reactions is that background contributions are strongly reduced com¬ 
pared to inclusive processes. A good example are searches for exclusive production of 
supersymmetric Higgs boson [3h5], anomalous boson couplings for 77 —> W + W~ 

[9] or for 77 —* 77 Mill. So far these processes are usually studied in the so-called 
equivalent-photon approximation (EPA) (for a description of the method, see e.g. M). 
Within the Standard Model the cross section for the pp —» ppW + W reaction is about 
100 fb at y/s = 14 TeV [13Q. Gluon-induced processes could also contribute to the exclu¬ 
sive production of W + W~ 1131 and W ± H T [14] via quark loops. 0 The corresponding 
cross sections are rather small mainly due to suppression of Sudakov form factors and 
the gap survival factor. The exclusive reactions could be also used in searches for neutral 
technipion in the diphoton final state [15] or dilaton [16]. Here a precise prediction of the 
cross sections is not possible as the model parameters are still unknown. 

Discovery of the heavy Higgs bosons of the Minimal Supersymmetric Standard Model 
(MSSM) [fl7[ -fT9ll or more generic Two-Higgs Doublet Models (2HDMs) (see e.g. [12(51. [ 2 lh 1 
poses a special challenge at future colliders. One of the international projects currently 
under consideration is the Future Circular Collider (FCC) ||22tl . The Higgs sector in 
both the MSSM and 2HDM contains five states: three neutral [two CP-even (h, H) and 
one CP-odd (A)] and two charged ( H + , H~) Higgs bosons. In general, either h or H 
could correspond to the SM Higgs. The charged Higgs boson pair production in the 
77 —* H + H~ mode was considered in |0-S]- In general, the higher-order corrections 
to the 77 —» H + H~ subprocess decrease the tree-level total cross section by about a few 
percent; see |26j,[27|]. Also the associated production 77 —> H ± W T was discussed in the 


literature [28]. For a more extensive discussion of charged Higgs boson production at the 


FHC and IFC, see 

There are also extensive phenomenological studies on charged Higgs boson(s) pro¬ 


duction at the FHC in the inclusive reactions via the partonic processes 
m H ± < m t — m^, the charged Higgs boson can be produced in t —* bH + and t 


31Q. If 
bH~ 


decays from the parent production channel pp —> tt, which would compete with the SM 
process tt —* bW + bW~. The dominant decay channels in this mass range are H± —> rv T 
and H ± —» cs(cs). In the case of a heavy charged Higgs with m H ± > mt — mj,, there are 
three major mechanisms: 

(a) Associated production with a top quark via the partonic processes qq,gg —» tbH ± 
[32-38] as well as through the gluon-bottom fusion gb —» tH ± |39 h 45(1 . The sequen¬ 
tial decay H + —> tb is known as a preferred channel. But signals in these processes 
appear together with large QCD backgrounds. The H ± —> W ± H/W ± A —> W^bb 
channels were analyzed in [}46ii47 ]. In the latter paper, the W ± tt decay channel was 
also considered. Recently, the H ± —* W ± (H 0 ^ S —> bb) decay channel for a SM-like 
Higgs was studied in 11481,14911 . This decay channel can be particularly important 
when charged Higgs is produced through the pp —> tH ± processes. 

1 An attractive channel is the associated production of a charged Higgs boson with a iW boson via 77 
fusion. Since there is no H + W r 7 couplings the 77 —» H + W 1 associated production process have no 
tree-level contribution in the 2HDM and in the MSSM and occurs only at one-loop level in the lowest 
order. 
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(b) Associated production with a boson through the qq,gg —» H ± W =F subprocesses 

|50l-l6Cll and associated production of a charged Higgs boson with a CP-odd Higgs 
boson, i.e. qq —* H ± A, was studied in 1 6ll. I62M . 


(c) Charged Higgs boson pair production via qq,gg —> H + H~ [63466], bb —* H + H~ 
[67]] subprocesses or in association with bottom quark pairs qq,gg —* bbH + H~ [68, 
l69[1 . For more recent studies, see 0:7l]]. 


The cross sections for the inclusive reactions strongly depend on the model parame¬ 
ters, such as tan /3 = vilv\, the ratio of the vacuum expectation values of the two Higgs 
doublets, and others. A program on how to limit the relevant parameters, based on the 
collider searches and data from B factories, was presented, e.g., in [721]. Another impor¬ 
tant ingredient of the model is the mass of the charged Higgs boson. In the MSSM the 
relation between the masses of the charged Higgs boson and CP-odd Higgs boson in 
lowest order is given by ra^ ± = m 2 A + mL ± 0 (for reviews and details, see, e.g. | 17Ul9 1). 

Several experimental searches already placed limitations on the mass of the charged 
Higgs bosons. There is a direct limit of m H ± > 78.6 GeV from the LEP searches f74fl 
by its decays H ± —> rv T and H ± —> cs(cs). At hadron colliders, the search procedures 
for a charged Higgs boson differ in term of its mass range. At the Tevatron the searches 
were mainly focused on the low mass range m H ± < mt which can put a constraint to the 
2HDM (as an example) on the small and large tan /3 regions for a charged Higgs boson 
mass up to 160 GeV |i75f] . Recent searches at the LHC [Za-|82] provide new limitations on 
the model parameters. However, still a possible span of parameters is rather large. For 
example, in the latest searches ATLAS and CMS put limits on the product of branching 
fractions BR(t —» H + b ) x BR(H + —» rv T ), but there are no model-independent limits on 
the H ± mass. The observed limits are reinterpreted in some MSSM scenarios, with mass 
limits around 140 — 160 GeV that depend somewhat on tan /L But in other models such 
as type-12HDM, the limit may be weaker. 

Other experimental bounds on the charged Higgs mass come from processes where 
the charged Higgs boson enters as a virtual particle, i.e. participates in loop diagrams. 
It is well known that in the type-II 2HDM, where the up- and down-type quarks and 
leptons couple to different doublets, the b —* S 7 transitions imposes a strong constraint 
on the Higgs boson mass m H ± A 300 GeV. In the type-I 2HDM, instead, all fermions 
couple to the same doublet and there is no such strong ^-physics constraint (the MSSM is 
also less sensitive to radiative corrections). The flavor constraints on the Higgs sector are, 
however, typically model dependent. A detailed analysis of precision and flavor bounds 
in the 2HDM can be found, e.g., in ll83ll . 

In the present analysis we wish to concentrate on exclusive production of charged 
Higgs bosons in proton-proton collisions proceeding through exchange of two photons. 
In Fig. |l] we show basic diagrams contributing to the pp —? ppH + H~ reaction. The 
coupling of photons to protons is usually parametrized with the help of proton electro¬ 
magnetic form factors: Gf (electric), Gm (magnetic) or equivalently Fi (Dirac), F 2 (Pauli). 
We wish to discuss the dependence on the form factors of several differential distribu¬ 
tions. In contrast to inclusive processes discussed above, the considered here exclusive 
reaction is free of the model parameter uncertainties, at least in the leading order, except 
of the mass of the charged Higgs bosons. 


2 This is particular to the MSSM in lowest order (is modified by one-loop radiative corrections l73ll l and 
does not hold in 2HDMs or in, e.g., the Next-to-Minimal Supersymmetric Standard Model (NMSSM). 
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FIG. 1. Born diagrams for exclusive production of pairs of charged scalar particles via photon- 
photon exchanges. 


Our paper is organized as follows. In Sec. HJ we discuss formalism of the pp —» 
ppH + H~ reaction both in the equivalent-photon approximation (EPA) in the momen¬ 
tum space commonly used in the literature and in exact 2 —* 4 kinematics. In Sec. HIT] we 
present numerical results for total and differential cross sections. In Sec. lIII Al we compare 
results obtained in the exact 2 —» 4 calculation and those obtained in EPA. We present not 
only estimation of the total cross section but also several differential distributions im¬ 
portant for planning potential future experimental searches. In addition, we discuss the 
role of absorption corrections commonly neglected for two-photon initiated processes. 
Finally, we also consider diffractive exclusive production of the H + H~ bosons through 
an intermediate recently discovered Higgs boson. Diffractive contribution is discussed 
in Sec. lIII Bl 


II. FORMALISM 

We shall study exclusive production of H + H~ in proton-proton collisions at high en¬ 
ergies 

p(Pa, A«) + p(p b ,h) -A p(pi, Ai) + H + (p 3 ) + FT(p 4 ) + p(p 2 ,A 2 ), (2.1) 

where p a ^, py 2 and A a ^, Ay 2 = denote the four-momenta and helicities of the pro¬ 
tons, and p 3/ 4 denote the four-momenta of the charged Higgs bosons, respectively. In the 
following we will calculate the contributions from the diagrams of Fig. 2] 


A. Equivalent-photon approximation 


Similar processes are treated usually in the equivalent-photon approximation (EPA) in 
the momentum space, see e.g. |13t 15ll . @ Only very few differential distributions can be 
obtained in the EPA approach. In this approximation, when neglecting photon transverse 
momenta, one can write the differential cross section as 


da 


dy 2 ,dy^d 2 p t H 16 n 2 s 2 


x 1 f(xi)x2f(x 2 )\M 7 ^ H + H - 


( 2 . 2 ) 


1 An impact parameter EPA was considered recently in [8; 
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where s = sx\X 2 and f(x)'s are an elastic fluxes of the equivalent photons (see e.g.[ 12 :]) as 
a function of longitudinal momentum fraction with respect to the parent proton defined 
by the kinematical variables of the charged Higgs bosons. 


X\ 


m t H 

V~s 


( e y 3+ gy4) / 


* 2 =^(e 3/3 + e XJ4 ), m tH = ^/\p tH \ 2 + m 2 H (2.3) 


with mtH being transverse mass of the H ± boson(s). Above \AA\ 2 is the 77 —» H + H 
amplitude squared averaged over the photon polarization states. 

The photon flux f(x) is given by the formula [12] 


f ( x ) = - , 
x Jq 2 


1 rQmax 0L em dQ 2 


n Q 2 


(1 - x) 1 - 


) D(Q 2 ) + yC(Q 2 ' 


(2.4) 


where the spacelike momentum transfer squared Q 2 = — q 2 = — t ^ 0 El and the photon 
minimal virtuality allowed by kinematics Q 2 njn = x 2 m 2 ( 1 — x)~ l . The coefficient func¬ 
tions C and D are determined by the electric and magnetic form factors of the proton: 

-1 


C(Q 2 ) = G z „m. D(Q Z ) = (4miG z E (Q z ) + Q z G z m (Q 


2 \ — 


ipC 2 ‘ 




Am 2 + Q 2 


, (2.5) 


where the Ge and Gm form factors are related to Dirac (Fi) and Pauli (F 2 ) form factors by 


G £ (Q 2 ) = fi(Q 2 ) - jLhQ 2 ), G m (Q 2 ) = Fi(Q 2 ) + F 2 (Q 2 ) 


( 2 . 6 ) 


Using the standard dipole parametrizations of the Sachs form factors (see, for instance, 
chapter 2 in I 86 lh 


G e (Q 2 ) = G d (Q 2 ) , G m (Q 2 ) = ^G d (Q 2 ) , 

V-N 

-2 
Q 2 

"‘D 


Gd(Q 2 ) = 1 + 


m 


, m 2 D = 0.71 GeV 2 , 


(2.7) 

( 2 . 8 ) 


D 


where Gd is the so-called dipole form factor, = 2.7928, }i p and pi ,v are the anomalous 
proton magnetic moment and the nuclear magneton, respectively, we obtain 

f '( Q2 >=( 1 + ^)( 1 + ^) Gd < q U <Z9) 

Gd(q2> ' (Z10) 

We shall use the parametrizations in the following analysis. 


4 Here we discuss the collinear EPA approach, that is, the photon transverse momenta qt 
proach including transverse momenta of photons was discussed recently in j&S]. 


0. An ap- 
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B. Exact kinematics 


In the present studies we perform, for the first time, exact calculations for the consid¬ 
ered exclusive 2 —» 4 process (12.11) . In general, the cross section can be written as 

j _ (^) 4 rrj- -[2 d 3 p 1 d 3 p 2 d 3 p 3 d 3 p A 

lL ' 2s 1 PP^PP H+ H-\ (2n 3 )2E 1 (2 tc 3 )2E 2 (2tt 3 )2E 3 (2tt 3 )2E 4 

xd 4 (E a + E b -p 1 -p 2 -p 3 -p 4 ) / (2-11) 


where energy and momentum conservations have been made explicit. The formula is 
written in the overall center-of-mass frame. Above |A1| 2 is the 2 —>■ 4 amplitude squared 
averaged over initial and summed over final proton polarization states. The kinematic 
variables for the reaction (12 .1 1) are 


S = (pa + Pb) 2 , S34 = M 2 H+H - = (p 3 + p 4 ) 2 , 

h = q\, h = q 2 ’ qi = Pa- p\, qi = Pb - pi • ( 2 . 12 ) 


Our calculations have been done using the VEGAS routine [87] and checked on an 
eight-dimensional grid H. The phase space integration variables are taken the same as 
in Ref.[88], except that proton transverse momenta p 3t and p 2t are replaced by = 
logio(pif / pot ) and £2 = logio(p2</po<)/ respectively, where = 1 GeV. The main ingredi¬ 
ents of the model are the amplitudes for the exclusive process. 

The Born amplitudes for the process (I2.1[) are calculated as 


M B Z% h A 2 mH-(tl,t2) = (2.13) 

where D /(y (f) = —igji V /1 is the photon propagator. Using the Gordon decomposition the 
7 pp vertex takes the form 


= eu(p', A') (V'fi (f) + ^—o-V v (p' - p) v F 2 (t)^j u(p. A) 

= eu(p',A') ^(Fi(f) + F 2 (t)) 7 f ‘ ~ 2 ^(P' + u(p,A) , ( 2 -14) 

where u(p, A) is a Dirac spinor and p, A and p’ , A/ are initial and final four-momenta and 
helicities of the protons, respectively. 

In the high-energy approximation, at not too large |f|,E|one gets the simple formula 

'm (p’ + py, (2.15) 

which is very convenient for the discussion of the proton spin-conserving and the pro¬ 
ton spin-flipping components separately. It is easy to see that in the approximation [see 
Eq. (12.151) 1 the cross section contains only terms proportional to Ff (t\)Fj (t 2 ) and no mixed 


5 The details on how to conveniently reduce the number of kinematic integration variables are discussed 
in H . 

6 We show how good the approximation is in Figs. 9,10,12,13. 
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terms proportional to F 2 (fi)Fi(f 2 )-P 2 (^ 2 )r etc. In exact calculations [with spinors of pro¬ 
tons, see Eq. (I2.14[) l. there is a small contribution of the mixed terms. This will be dis¬ 
cussed when presenting our results. 

The tensorial vertex in Eq. (I2.13|) for the 77 —> H + H~ subprocess is a sum of three- 
level amplitudes corresponding to t, u and contact diagrams of Fig. [T[ respectively, 

V v 7 ^ h+h - = Vt lV2 + Vu lV2 + vr 2 


= le 


pf — mf 


+ie 


Vl ~ m H 


■{qi - P4 + P3) Vl (qi - 2p 4 ) V2 

(qi - 2p 4 ) Vl (^i - p 4 + P3) V2 - 2 ie 2 g VlVl , 


(2.16) 


where p] = (q 2 ~ p 4 ) 2 = {qi ~ P3 ) 2 and p 2 = - p 4 ) 2 = (q 2 - p 3 ) 2 - There are strong 

cancellations between the three contributions. 

A complete calculation for exclusive H + H~ production in pp collisions, in addition 
to the 77 exchange, must take into account more diagrams than those of Fig. |T[ We can 
have the 7 Z, Z 7 , and ZZ exchanges. The corresponding amplitudes can be obtained by 
substitution of the photon propagator and the 7 pp vertex [see (12.131) and (12. 14[) 1 by the Z 
boson propagator and the Zpp vertex [89], 


rifZfAo = A-s(pU') (y'fT c (t) + + 7 " 75 G" c (t)) u(p, A), 


SwCw 


(2.17) 


where we use the shorthand notation cw = cos % = sin 9yj, 6w is the Weinberg 
mixing angle. The 7 H + H~ and 77 H + H~ coupling constants in (12.161) read: 


gZH+H- 

gryZH+H- 

gZZH+H- 


ie 


2cwSw 


( C W “ s w) / 


ze 


CwSw 

„2 


( C W S^y) , 


ze 


2 \2 


2 c 2 s 2 


( c w _ s w) 


(2.18) 


The neutral current form factors appearing in (12.171) related to the vector part can be 
related to electromagnetic form factors [see (12.91) and (12.101) . Fi >2 (t) = Fj’ 2 (f)], 


Fl N 2 C (t) = \ (r£ 2 (f) - 72 ( 0 ) - 28^(0 - lf{, 2 (f), (2.19) 

where |F”(f)| <C Fj’(f) for small \t\. The form factor related to axial-vector neutral 
current is related to the familiar charge current axial-vector form factor: 


G« c (f) = lc^(t) - |G^(f). (2.20) 

Since the strangeness form factors F| 2 arid G S A are poorly known and small in the follow¬ 
ing we shall neglect them. 
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C. Absorption corrections 


The absorptive corrections to the Born amplitude (12.131) are added to give the full phys¬ 
ical amplitude for the pp —» ppH + H~ reaction: 

M pp ^ ppH+H - = M®°™ ppH+H - + • ( 2 - 21 ) 

Here (and above) we have for simplicity omitted the dependence of the amplitude on 
kinematic variables. 

The amplitude including pp-rescattering corrections between the initial- and final- 
state protons in the four-body reaction discussed here can be written as 

M 2a^aXh* h (S'Pv'Pit) =^r s J ^M KX ^ KK {s,-k\) 

X -^A°A(^AiA 2 H+H-( S '^lf' Pit ) ' 

( 2 . 22 ) 

where p lt = p lt — kt and p 2t = p 2t + kt■ Here, in the overall center-of-mass system, p lt 
and p 2t are the transverse components of the momenta of the final-state protons and kt is 
the transverse momentum carried by additional pomeron exchange. M. pp ^ pp (s, —kf ) is 

r\ 

the elastic pp-scattering amplitude for large s and with the momentum transfer t = —kj. 
Here we assume s-channel helicity conservation and the exponential functional form of 
form factors in the pomeron-proton-proton vertices. 

We shall show results in the Born approximation as well as include the absorption 
corrections on the amplitude level. This allows us to study the absorption effects dif¬ 
ferentially in any kinematical variable chosen, which has, so far, never been done for 
two-photon induced (sub)processes. 


III. RESULTS 


A. Electromagnetic process 


In this section we shall present results of our calculations for the pp —» ppH + H~ 
reaction (12 . T |) calculating from the diagrams of Fig. [I] Let us start our presentation by 
presenting the total cross section for y/s = 14 TeV (LHC) and y/s = 100 TeV (FCC) and 
for various charged Higgs mass values. In Table U we show cross sections in fb without 
and with (results in the parentheses) the pp-rescattering corrections. The smaller the 
values of m H ±, the larger are those of cross section 0. The values of the gap survival 
factor (S 2 ) for different masses of H 1 * 1 bosons m H ± = 150,300,500 GeV are, respectively, 
0.77,0.67,0.57 for y/s = U TeV (LHC) and 0.89,0.86,0.82 for y/s = 100 TeV (FCC). In 
contrast to diffractive processes, the larger the collision energy, the smaller the effect of 
absorption. We have checked numerically that the cross section contributions with the 
7 Z, Z 7 , and ZZ exchanges are very small compared to the 77 contribution and will be 
not presented explicitly in this paper. 


7 We wish to note on the margin that the cross section for pair production for doubly charged (Higgs) 
bosons, e.g. H ++ H , would be 16 times larger 03 in the leading-order approximation considered 
here. The doubly charged Higgs bosons are expected in models that contain a Higgs boson triplet field. 
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m H ± (GeV) 

150 

300 

500 

Clhc (fb) 
CTcc (fb) 

0.1474 (0.1132) 
1.0350 (0.9236) 

0.0119 (0.0080) 
0.1470 (0.1258) 

0.0014 (0.0008) 
0.0303 (0.0249) 


TABLE I. Cross sections in fb for the pp —> ppH 1 H reaction through photon-photon exchanges 
without and with (results in the parentheses) the absorption corrections for two center-of-mass 
energies y/s = 14 TeV (LHC) and y/s = 100 TeV (FCC) and various charged Higgs bosons mass 
values. The calculations was performed for exact 2 —> 4 kinematics and with the amplitudes in 
the high-energy approximation, see Eq. (12.15b . 



FIG. 2. Distribution in the auxiliary variables </\ or £2 at y/s = 14 TeV (left panel) and 100 TeV 
(right panel). The online red solid lines represent the calculation of exact amplitude (including 
spinors of protons). The black upper and lower long-dashed lines correspond to calculations in 
the high-energy approximation (12.151) without and with the absorption corrections, respectively. 


In Fig.|2]we show a distribution in an auxiliary integration variable(s) £ 1/2 = log 10 (pi/ 2 f/I GeV). 
If protons are measured, the distributions in Fig. |2|can be measured too. Here and in the 
following, we discuss the differential distributions for one selected mass of H ± . For ex¬ 
ample, we shall assume m H ± = 150 GeV, which is rather a lower limit for the charged 
Higgs bosons. The general features of the differential distribution for heavier masses are, 
however, similar. We compare results without (the upper long-dashed lines) and with 
(the lower long-dashed lines) absorption corrections due to the pp interactions. 

The rapidity distribution for the charged Higgs bosons is shown in Fig. [3j The larger 
center-of-mass energy the broader the rapidity distributions. 

In Fig. H] we show invariant mass distribution of the H + H~ subsystem in a broad 
range of the invariant masses. We compare results for the exact kinematics and for the 
EPA calculations. Please note that for the EPA, the invariant mass of the diHiggs system 
is given by M H + H - ~ sx\X 2 - 

In Fig. [5] we show decomposition into helicity components of the cross section in the 
two-Higgs invariant mass and in the rapidity of one of the charged Higgs bosons. Here 


9 














FIG. 3. Rapidity distribution of charged (Higgs) bosons at ,/s = 14 TeV (left panel) and 100 TeV 
(right panel). The meaning of the lines is the same as in Fig. [2] The short-dashed (online green) 
lines represent results of EPA. 
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FIG. 4. DiHiggs boson invariant mass distributions at \/s = 14 TeV (left panel) and 100 TeV (right 
panel). The meaning of the lines is the same as in Fig.|2| The short-dashed (online green) lines 
represent results of EPA. 


we use the formula (12.151) for the 7 pp vertex which is very convenient for the discussion of 
the proton spin-conserving (the Dirac form factor (12.91) only) and the proton spin-flipping 
(the Pauli form factor (12.10)) only) components separately. 

In Fig. [ 6 ] we show the dependence of absorption on M H + H -. This is quantified by the 
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FIG. 5. DiHiggs boson invariant mass (left panel) and Higgs boson rapidity (right panel) distri¬ 
butions for v/i = 14 TeV in the Born approximation and for the amplitudes given by Eq. (12.15b . 
The double spin-conserving contribution (cc) is show by the long-dashed line, the double spin¬ 
flipping contribution (ff) by the dotted line, and the mixed contributions (cf) and (fc) by the dashed 
lines. The solid line represents the sum of all the contributions and corresponds to the upper long- 
dashed lines in Figs. |3] and [4] (left panels). 


ratio of full (with the absorption corrections) and Born differential cross sections 


(S 2 (M h+h -)) 


(icjBorn + absorption / dM H + H 

d<r Born /dM H+H - 


(3.1) 


The absorption effects due to the pp interaction lead to large damping of the cross section 
at the LHC and relatively small reduction of the cross section at the FCC. This result must 
be contrasted with typical diffractive exclusive processes where the role of absorption 
effects gradually increases with the collision energy. 

In Fig. [7] we show the ratio of the cross section for all (Fi, F 2 ) terms included in the 
amplitude to that for F\ terms only both for the exact 2 —» 4 kinematics and for the EPA 
calculations. Here for consistency we have neglected the interference effect between the 
electromagnetic form factors in the EPA approach. At large invariant masses of Mp j+h- 
the ratio for exact calculation is much smaller than that for EPA. This suggests that EPA 
overestimates the spin-flipping contributions. 

Let us discuss now a subtle effect of the interference of terms proportional to F\ and 
F 2 . To quantify the effect, let us define the following quantities: 


da mcoh = ^(FpFO + da(F 1 ;F 2 ) + dcr(F 2 ;F 1 ) + da(F 2 ;F 2 ) , (3.2) 

dcr C0h = da(F ll F 2 ;F 1/ F 2 ), (3.3) 


where dcr(Fp,Fj ) means the cross section when, at one proton line, only the F, term is 
taken into account and, at the second proton line, only the Fy term is taken into account. 
dcr coh represents the cross section where all terms are coherently included. In Fig. [ 8 ] we 
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FIG. 6. The dependence of the gap survival factor due to pp interactions on M H + H - for exact 
2 —> 4 kinematics at = 14 TeV (left panel) and 100 TeV (right panel). This is quantified by the 
ratio of full (including absorption) and Born differential cross sections ( S 2 (Mh+h -)) ( 13.11 1. 


T3 


Ph 

Id 

T3 



M h+h (GeV) 


M h+h (GeV) 


FIG. 7. Ratio of the cross section when all (F\, F 2 ) terms are included to that with Fi terms only 
for exact 2 —> 4 kinematics (black long-dashed lines) and EPA (online green short-dashed lines) 
calculations at \fs = 14 TeV (left panel) and >Js = 100 TeV (right panel). The lowest long-dashed 
lines represent results with the absorption corrections. 


show the relative corrections (( da coh — dcr incoh ) /dcr coh ) coming from the interference effect 
between different terms in the amplitude. We see from Fig. [8] that for ^ = 2 TeV 

the total cross section from the calculation using exact amplitude (including spinors of 
protons) is modified by ~ 10% at \/s = 14 TeV, while at y/s = 100 TeV only by ~ 1%. The 
smallness of the effect causes the effect of the fluctuations in our Monte Carlo approach. 
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FIG. 8. The relative corrections in percent as a function of M H + H - for exact Born amplitude in¬ 
cluding spinors of protons (online red solid lines) and for EPA (online green short-dashed lines) 
calculations at ^/s = 14 TeV (left panel) and 100 TeV (right panel). 


The relative corrections for the EPA approach are somewhat larger. 

In Fig. [9] we present distributions in charged Higgs boson transverse momentum pt,H> 
i.e., pt r H+ or Pi, 11 ■ While at low (Higgs) boson transverse momenta the EPA result is very 
similar to our exact result for all spin components, some deviations can be observed at 
larger transverse momenta. This is consistent with the similar comparison for the distri¬ 
butions in invariant mass (for the process under consideration large transverse momenta 
are related to large invariant masses). 

If forward/backward protons are measured, then distributions in four-momentum 
transfers squared (t = t\ or £ 2 ) can be obtained and relevant distributions shown in Fig.HOl 
can be constructed. The absorption effects due to the pp interactions are stronger for large 
values of \t\. 

In Fig. [IT] we show a decomposition of the cross section into helicity components as 
a function of momentum transfer(s) squared. The proton spin-conserving contribution 
related to the Dirac form factor(s) clearly dominates at very small |fi| or 1 12 1 - At larger \t\ 
the proton spin-flipping contribution related to the Pauli form factor(s) becomes impor¬ 
tant as well. The double spin-flipping contribution (ft) vanish at |ti| = | f 2 1 =0, while the 
mixed contributions (fc) and (cf) vanish at 1 1 11 =0 and | f 21 = 0, respectively. 

Fet us consider now azimuthal correlations between outgoing particles. In Fig. [12] 
we show correlations between outgoing protons. We emphasize the dip at (p pp — zr /2 
which is a consequence of the couplings involved in calculating the 77 —> H + H~ matrix 
element(s). 

The correlation between outgoing Higgs bosons is shown in Fig. [13[ The bosons are 
produced preferentially back-to-back which can be understood given small transverse 
momenta of virtual photons compared to transverse momenta of the Higgs bosons. 
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FIG. 9. The distributions in the Higgs boson transverse momentum at ^/s = 14 TeV (left panel) 
and 100 TeV (right panel). The meaning of the lines is the same as in Fig. |21 The online green 
short-dashed lines represent results of EPA. 



FIG. 10. Distribution in momentum transfer(s) squared (h or ti ) at = 14 TeV (left panel) and 
100 TeV (right panel). The meaning of the lines is the same as in Fig.|2| 


B. Diffractive process 


So far we have considered a purely electromagnetic process, the contribution of which 
is model independent. The corresponding cross section turned out to be rather low. 
Therefore, one could worry whether other processes might not give a sizeable contribu¬ 
tion, comparable to the photon-photon exchanges. One such candidate is the diffractive 
mechanism discussed, e.g., in the context of exclusive Higgs boson production [93—96]. 
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FIG. 11. Distributions in momentum transfer squared t\ (left panel) and tj (right panel) at 
\fs = 14 TeV in the Born approximation and for the amplitudes given by Eq. (12.15b . The dou¬ 
ble spin-conserving contribution (cc) is show by the long-dashed line, the double spin-flipping 
contribution (ff) by the dotted line, and the mixed contributions (cf) and (fc) by the dashed and 
dot-dashed line, respectively. The solid line represents the sum of all the contributions and corre¬ 
sponds to the upper long-dashed line in Fig. [lOl 


In the present case the mechanism shown in Fig. [14] seems an important candidate. The 
g*g* —> H + H~ hard subprocess amplitude through the f-loop and s-channel SM Higgs 
boson (h°) is given by 


V g . g ^m„- - V gi ^„ - - m ,' + — gi /H + H - 


(3.4) 


and enters into M. 


pp^>ppH+H~ 


invariant 2 —> 4 amplitude for the diffractive process as in 


[hiH]. The triple-Higgs coupling constant ghn+H- is, of course, model dependent. In the 
MSSM model it depends only on the parameters a. and p. In the general 2HDM it depends 
also on other parameters such as the Higgs potential A-parameters or masses of Higgs 
bosons. How the coupling constant depends on parameters of 2HDM was discussed, 
e.g., in 197- lodl . In Fig[T5]we show as an example the coupling as a function of tan p and 
a. — p for MSSM (left panel) and 2HDM (right panel). In the latter case we have used 
a relation given in Ref. |J99|] while the formula for the MSSM can be found e.g. [17]. The 
ghH+H- coupling constant in the MSSM case does not exceed 50 GeV (to be compared e.g. 
t° ghhh ~ 194 GeV in the Standard Model). The coupling constant in the case of 2HDM 
can be, in general, very large. Recent data obtained at the THC in the last three years 
put stringent constraints on a. and p as well as on masses of the, thus far, unobserved 
Higgs bosons (some examples of such analyses can be found in |83, 10(1 101 [. The THC 
experimental data allow for two regions in the [tan p, (ft — a)] plane [83,100, 10211 . One of 
them p — a. ~ 7 t/2 is the so-called "alignment limit". The second one is more difficult to 
characterize. In the present analysis we focus on the alignment region which means the 
lightest CP-even Higgs h is what has been found at the THC with m /, ~ 125 GeV [1031] . 
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FIG. 12. Distribution in relative azimuthal angle between outgoing protons at >/s = 14 TeV (left 
panel) and 100 TeV (right panel). The meaning of the lines is the same as in Fig. [2] 




FIG. 13. Distribution in relative azimuthal angle between outgoing charged (Higgs) bosons at 
\/s = 14 TeV (left panel) and 100 TeV (right panel). The meaning of the lines is the same as in 

Fig. El 


Experimental data allow for some deviations from the ft — ol = n/ 2. As can be seen from 
Fig. [15] a small deviation from this limit can modify the coupling constant considerably. 
The analysis in [83] suggests that m H ± ~ tha and we keep such a relation throughout 
our analysis. A deviation from such a relation would increase the discussed coupling 
constant. 

In Fig. [16] we show the dependence of the coupling constant on masses of charged and 
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FIG. 14. The diffractive mechanism of the exclusive charged Higgs bosons production through 
the intermediate CP-even neutral recently discovered Higgs boson. The absorption corrections 
due to pp interactions (indicated by the blob) are relevant at high energies. 
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FIG. 15. The triple-Higgs coupling constant \g hH + H - \ in the [tan /3, (/3 — a)] plane for the MSSM 
(left panel) and the type-II 2HDM with m H ± = in a = 300 GeV (right panel). Here we consider 
the type-II 2HDM in which the scalar potential parameters are Re( A 5 ) (explicit Z 2 symmetry) and 

3-6 = A7 = 0. 


CP-odd Higgses within 2HDM. @ A minimal value appears when ~ m H ±. When we 
relax this condition the coupling can be even as large as 1000 GeV. This is consistent with 
the limits of the allowed region in [ 100h . Summarizing, ghH+H- ' n the 2HDM is limited 
to 64 GeV < gj,i-{ 1 a < 1000 GeV. The corresponding couplings in the MSSM are smaller 
than 50 GeV.@ 

In Fig. [17] we show corresponding results for the diffractive contribution for \/s = 14 
TeV (left panel) and sjs = 100 TeV (right panel) both the lower and upper limit of the 
2HDM triple-Higgs coupling for m H ± = 150 GeV. In the calculation we have included 
the "effective" gap survival factor (S 2 ) = 0.03 typical for the considered range of ener¬ 
gies. The cross section for the exclusive diffractive process is much smaller than that for 


8 We emphasise again that in the MSSM in the lowest order we have ?«T ± = m 2 A + m 2 ±. 

9 It has been shown, e.g., in li57ll that in some regions of the parameter space of 2HDMs the associated 
production cross section can be enhanced compared with the MSSM by orders of magnitude. This is a 
similar process to that discussed in our paper. 
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FIG. 16. | ghH+H~ I f° r the 2HDM as a function of m H ± and for the alignment limit f J > — a = n/2. 




FIG. 17. DiHiggs boson invariant mass distributions at n/s = 14 TeV (left panel) and 100 TeV 
(right panel). The upper lines represent the 77 contribution. We also show contribution of the 
diffractive mechanism (the shaded area) for the MSTW08 NLO collinear gluon distribution [104] 
and ghH+H- = 100 (1000) GeV for the lower (upper) limit. 


77 mechanism both for LHC and FCC. In addition to the result for the 2HDM set of pa¬ 
rameters (alignment limit), we also show result with the upper limit ghH+H- = 1000 GeV. 
With such a big coupling constant, the contribution with the intermediate neutral Higgs 
boson h° dominates over the contribution of boxes for tan < 20. Therefore, the upper 
limit also effectively includes the box contributions discussed in the context of inclusive 
PP “V (88 “V H + H ~) processes [63-65]. 
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C. A comment on possible experimental studies 


So far we have calculated the cross section for the pp —» ppH + H~ reaction. If one 
wishes to identify the reaction experimentally, one should measure the decay products 
of H± bosons. The branching fractions to different channels depend on the model pa¬ 
rameters ( m H ± , tan ft, etc.). For low masses of H ± (m H ± < 150 GeV) it is expected that 
the rv T and cs are the dominant channels. In the case of heavy charged Higgs (with 
m H ± > 200 GeV) the tb(tb) or W ± h channels are expected to be the relevant ones. 

In the first case (light H ± ) t + t~ could be measured in addition to the forward/backward 
protons. The emission of neutrinos leads to a strong imbalance between proton-proton 
missing mass and M T + T -. This should help to eliminate the pp —» ppr + T~ reaction, 
but this requires dedicated Monte Carlo studies, including perhaps the pp —» ppr + r~j 
process. The pp —» ppW + W~ and pp —» ppH^W 11 reactions may lead to a similar final 
state. Although the branching fraction W + —* r + v T or W“ —» r~v T is only about it 
is expected to be a difficult irreducible background because of the relatively large cross 
section for the pp —* ppW + W~. In principle, the cscs (four jets) final channel is also 
attractive as, in this case, one may check extra conditions M^i — m w + >10 — 20 GeV and 
Mqq' — m w - >10 — 20 GeV to exclude the W + W~ contribution. The pp —> ppH + W~ 
and pp —» ppH~W + processes may lead to similar final channels (r + v T r~v T or cscs ) but 
the corresponding cross sections are expected to be smaller (higher-order processes with 
loops). Mixed (leptonic + quarkish) final states could also be considered. 

In the second case (heavy H ± ), in general, both the t quark and b jet can be mea¬ 
sured. In contrast to the previous case we do not know about any sizeable irreducible 
background. But then the cross sections are rather small as discussed in the previous sec¬ 
tions. In the case of the H + H~ —» W + hW~h decay channel the actually measured final 
state can be rather complicated (e.g., qq'bb ifcjbb). Therefore, with experimentally limited 
geometrical acceptance it may be rather difficult to reconstruct the charged Higgs bosons. 

A detailed analysis of any of the final states considered here requires separate Monte 
Carlo studies including experimental geometrical acceptances relevant for a given exper¬ 
iment. This clearly goes beyond the scope of the present paper which aims to attract 
attention to potentially interesting exclusive processes. The Monte Carlo studies could 
be done only in close collaboration with relevant experimental groups. 


IV. CONCLUSIONS 

In the present paper we have studied in detail the exclusive production of heavy scalar, 
weakly interacting, charged bosons in proton-proton collisions at the LHC and FCC. In 
contrast to EPA our exact treatment of the four-body pp —* ppH + H~ reaction allows us 
to calculate any single particle or correlation distribution. 

Results of our exact (2—^4 kinematics) calculations have been compared with those 
for the equivalent-photon approximation for observables accessible in EPA. Rather good 
agreement has been achieved in those cases. However, we wish to emphasize that some 
correlation observables in EPA are not realistic, or even not accessible, to mention here 
only correlations in azimuthal angle between the outgoing protons or the charged Higgs 
bosons. We have predicted an interesting minimum at (ppp = 90° which is a consequence 
of the field theoretical couplings involved in the considered reaction. 
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We have analyzed in detail the role of the Dirac and Pauli form factors. In contrast to 
light particle production, the Pauli form factor plays an important role especially at large 
M h + h -, and related terms in the amplitude cannot be neglected. We see that the double 
spin preserving contributions are almost identical in both exact and EPA calculations 
(within 1%), but the spin-flipping contributions are in our calculation somewhat smaller. 

In the present paper we have studied, for the first time for the considered two-photon- 
induced reaction, the absorption effects due to proton-proton (both initial and final state) 
nonperturbative interactions. Any extra interaction may, at the high energies, lead to 
a production of extra particles destroying exclusivity of the considered reaction. The 
absorptive effects lead to a reduction of the cross section. The reduction depends on 
kinematical variables. A good example are distributions in four-momentum transfers 
squared. At small |fi | and \ t 2 \, the absorption is weak and increases when they grow. We 
have also found interesting dependence of the absorption on 

The relative effect of absorption is growing with growing A similar tendency 

has been predicted recently for the pp —* ppW + W~ in the impact parameter approach 
[84]. The impact parameter approach is, however, not useful for many observables stud¬ 
ied here. We have predicted that the absorption effects for our two-photon-induced pro¬ 
cess become weaker at larger collision energy which is in contrast to the typical situation 
for diffractive exclusive processes. Our study shows that an assumption of no absorption 
or constant (small) absorption effects, often assumed in the literature for photon-photon- 
induced processes, is rather incorrect and corresponding results should be corrected. 

In addition to calculating differential distributions corresponding to the 77 mecha¬ 
nism we have performed first calculations of the H + H~ invariant mass for the diffrac¬ 
tive KMR mechanism. We have tried to estimate limits on the ghn+H- coupling constant 
within 2HDM based on recent analyses related to the Higgs boson discovery. The diffrac¬ 
tive contribution, even with the overestimated | ghn+H-1 coupling constant, gives a much 
smaller cross section than the 77 mechanism. We have also made an estimate of the 
contributions related to 7 Z, Z 7 , and ZZ exchanges and found that their contributions 
are completely negligible. This shows that the inclusion of the 77 mechanism should be 
sufficient, and the corresponding cross sections should be reliable. 

Whether the pp — > ppH + H~ reaction can be identified at the LHC (run 2) or FCC re¬ 
quires further studies including simulations of the H ± decays. Two H± decay channels 
seem to be worth studying in the case of light H ± : H± —* t + v t ( t “ v t ) or H ± —» cs(cs). 
The first decay channel may be difficult due to a competition of the pp —> ppW + W~ re¬ 
action which can also contribute to the t + t “ channels. The combined branching fraction 
is about 0.11 2 = 0.0121 (two independent decays) which is not so small given the fact that 
the cross section for the W + W“ production is much bigger than that for H + H~ produc¬ 
tion. In the second case (four quark jets), one could measure invariant masses of all dijet 
systems to reduce the W + W~ background. In the case of the heavy H ± Higgs boson, the 
H 1 * 1 —> tb{tb) decay can be considered. In principle, both the t quark and b jet can be 
measured. In this case we do not know about any sizeable irreducible background. 

The reaction considered in this paper is a prototype for any two-photon-induced pro¬ 
cess. In the future we wish to also consider the pp —> ppW + W~ reaction where similar ef¬ 
fects may occur. This reaction was proposed to search for the anomalous triple or quartic 
boson coupling. Effects beyond the Standard Model are expected at rather large invari¬ 
ant masses M w + w -, where we have found strong absorptive corrections. This exclusive 
reaction is, however, more complicated due to the more complex couplings and spins in- 
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volved and due to weak decays of the two W bosons where strong spin-spin correlation 
effects are expected. 
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